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MULTIGRADED POINCARE SERIES FOR MIXED STATES OF 

TWO QUBITS AND THE BOUNDARY OF THE SET OF 

SEPARABLE STATES 

O ' DRAGOMIR Z. BOKOVIC 

O 

(N 

, Abstract. Let A^ be the set of mixed states and S the set of separable states 

*^^i of the two-qubit system. Its Hilbert space is the tensor product H = C Cgi C^, 

and the group of local unitary transformations is G = SU(2) X SU(2) (ignoring 
the overall phase factor). Let V be the algebra of real polynomial functions 
on the space of all hermitian operators of trace 1 on "H. Let V^ C P be the 
OO ■ subalgebra of G-invariants. We compute its multigraded Poincare series and 

verify that it is consistent with Makhlin's list of 18 invariants. By using the 
Cn ' recent result of Augusiak et al. we describe the boundary of S and show that 

^ , its intersection with the (relative) interior of At is a smooth manifold. 

o 

O ■ 1. Introduction 

f^ ' Let C" be the Hilbert space (of column vectors) with the inner product (xi ,X2) — 

x\x2 and let M„ be the algebra of complex n x n matrices, li X,Y £ Mn then 
Cini their inner product is defined by (Xi,X2) = tr(X|X2). Let su{n) denote the Lie 

algebra of SU(n), it consists of all traceless skew- hermitian matrices in M„. 
G I The quantum system that we consider is bipartite. Both parties are qubits, C^. 

It will be convenient to identify C^ with the tensor product Ti. = C^ ®c C^, and 
qh; M4 with M2 00 M2. For Xi,X2,Yi,Y2 G M2, the inner product of Xi (g) Yi and 

X2 Y2 is given by 



X 



{Xi<E>Yi,X2<^Y2) = {Xi,X2) ■ {Yi,Y2). 



H \ It follows that {X ® Yy — X'^ ® Y'^ and, in particular, the tensor product of 

hermitian matrices is again a hermitian matrix. 

The space of hermitian matrices H^ C M^ is the direct sum of the \D real space 
spanned by the identity matrix, /4, and the space H^^q = i5u(4). The space H^fi 
is the direct sum of three real subspaces: Vi = 7^2,0 ® I^i V2 = I2 ® H2fi, and 
V3 — i?2,o ® ^2,0- Any mixed state, p, of our quantum system can be written 
uniquely as the sum of four components; 

(1.1) p^-h+X®l2+l2(SY + Z, 

where X,Y e H2,o and Z eVs. 

Denote by G the group of local unitary transformations, SU(2) x SU(2), where 
we ignore the overall phase factor. Note that G acts on M4 in the usual manner: 
{g, Z) — !■ gZg^ and stabilizes the real subspaces H4 and H4^o. Moreover, each of 
the subspaces Vi, V2, V3 is a simple G-module. 
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2 D.Z. BOKOVIC 

Let V denote the algebra of real valued polynomial functions on H^.t) and V'^ 
the subalgebra of G-invariant functions. Note that 7"^ inherits the Z-gradation 
from V ■ The homogeneous polynomials of V'^ may be used to construct measures 
of entanglement of our quantum system. 

The knowledge of 7"-^ is important because of the following well known fact: 
Two states, say pi and p2, belong to different G-orbits iff there exists / S V'^ such 
that /(/Oi) ^ f{p2)- If this holds true for a subalgebra A C P'^, then we say that 
A is complete. 

The Poincare series (also known as the Hilbert series) of V'^ was computed by 
M. Grassl et al. |4i Section VI]. We have 

oo 

P{z) = Y,d\u,{Vf)z\ 
d=a 

where V^ is the space of homogeneous polynomial G-invariants of degree d. This 
Poincare series is a rational function of the variable z: 

(1.2) P{z) — o « T Q • 

(1-Z)'(l + Zf (1 + Z2)2(1+Z + Z2)3 

We remark that in T the formula contains an extra factor 1 — z in the denominator 
because the authors work with the space H^ instead of H/^^q- 
The Taylor expansion begins with 

P(Z) : 



Subsequently, Y. Makhlin jSi gave a simple construction of 18 invariants which 
generate a complete subalgebra of V'^ . This subalgebra is proper and he gives 
two additional invariants. However, it is still not established whether this enlarged 
subalgebra is in fact the whole algebra V'^ . (According to M. Grassl [H] this can 
be proved by using the relations for the invariant tensors of SU(2).) 

As a complement to the two papers just mentioned, we compute in Section|2]the 
Z'^-graded Poincare series of V . This multigraded series carries more information 
about the invariants than the simply graded one. 

Let us recall that a mixed state is called separable if it can be represented as a 
convex combination of product states, and otherwise it is called entangled. Denote 
by Ai the set of all states and by S the set of all separable states. It is well known 
that these two sets are compact and convex, and have nonempty interior. By dA4 
resp. dS we denote the boundary of Ai and S, respectively. 

In Section|3we show that there is a natural way of decomposing the boundary dS 
into two pieces. Each of these pieces is a portion of a real algebraic hypersurface 
in the ambient afline space. In particular, it follows from our results that the 
intersection of dS with the relative interior of A^ is a smooth manifold. 

2. The multigraded Poincare series 

The direct decomposition H^fi = Vi V2 ® V3 induces a Z'^-gradation on V which 
is preserved by the action of G. If the space of homogeneous invariants of degrees 
di, d2, ds (with respect to the coordinates of the three subspaces) has dimension 



1 + 3 


z^ + 2z^ + 10z'^ + 7z^ 


+ 29zS + 25z^ 


+ 73z^ + 74:z^ + 172z'^° 


+187 


z" + 381 2^2 + 431 2^3 


+ 785z" + 920 


z^^ + 1539 z^^ + 1827 z" 


+287 


8 2^^ + 3441 2^9 + 5151 


z20 + 6185z2i4 


8887 z22 + 10666 z23 + ... 
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c, then this fact is recorded in the Taylor expansion of the multigraded Poincare 
series by the term ct'^H'^^tf'. 

Theorem 2.1. The multigraded Poincare series P(ti,t2,t3) of the algebra V of 
local unitary polynomial invariants of mixed states of two qubits is the rational 
function, whose numerator N{ti,t2,t:j) and denominator D{ti,t2,t3) are given by: 

N ^ 1- tits^ - t2t3^ + iitats^ + tit2t3^ + tiH2t3^ + tiia^^s^ + hhs^ + iitats* 

+t2 ts — ti ^2^3 — ti t2 ts — tit2 ts — ti ^2^3 — ^1^2 ^3 
— ti t2 ^3 — ti t2 t^ + ti t2 ^3 + ti t2 is — tl t2 t^ , 
D = (l-tl2)(l-t2')(l-t3')(l-ili2i3)(l-tli3') (1-^2*3') (1-^3') 
{l-t,W){l-t2W){l-t,^). 

The theorem is proved by using the weU-known Mohen-Weyl formula. Following 
the recipe from '^ , we obtain that 

, , ^ f f f I , da; dy dz 

P{.ti,t2,t3) = / / / ip(X,y,Z,ti,t2,t3) , 

(^'^^) J\z\=iJ\y\=iJ\x\=i X y z 

where 

^{x, y, z, t,,t2M) = ^, V zt.t2 t.) 

and 
V' = {l-ti)(l-t2f{l-t:if{l-tix){l-tix-'){l-t3xf{l-t3X-^f 

(1 - t3xy){l ~ t3xz){l - t3xyz){l - t3xy^'^){l - t3xz^'^){l - t^xy'^z'^) 

(1 - t3X-^y){l - t3X-^z){l - tsx-^yz) 

(1 - t3X~^y~^){l - t3.x-iz-i)(l - tsx-^y-^z-^) 

(1 - t2y)(l - t2yz){l - t2y-'){l - t2y-^z-^) 

(1 - t3y)(l - hyz){l - t3y-')(l - t:iy-^z-^) 

(1 - t2z)(l - t2z"')(l - t3^)(l - t:iz-^). 

The three integrations are to be performed over the unit circle in the counterclock- 
wise direction assuming that the \ti\ < 1. The computation was carried out by 
using Maple ^. 

After setting ii = t2 = ts = z in P{ti,t2,t^), the numerator and denominator 
acquire the common factor (1 + z^)(l — z^). After cancellation, we obtain exactly 
the Eq. p.2|l . i.e., we have P{z, z, z) = P{z). 

By expanding P(ti,t2,i3) in the Taylor series, we find that 

P = l + ti^+t2'^+t3'^+t3^+tit2t3 + 2t3-^ + 2t2%^+t2* + tit2t3^ 

+ 2ti^t-i^ + tiH2^ + ti^ + ta^ + t2^<3''^ +2tit2t3^ + tif2^t3 + h^ti^ 
+tih2t3 + 3 t3^ + 4 t2^t3^ + 2 t2%^ +t2^ + 2 tit2t3^ + tit2%^ + ^1^2^*3^ 
+4*12^3* + tiH2t3^ + Ati'^t2^t3'^ + tiH2^ + tih2t3^ + 2 tiHs^ 

To make the connection with the notation in Makhlin's paper, we may assume 
that the components of his vectors s and p are the linear coordinates on Vi and 
V2, respectively, and the nine entries [3ij of his matrix (3 are the coordinates on 
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V3. Then it is easy to check that his hst of invariants agrees with the information 
provided by the coefficients in the above Taylor expansion. For instance, the cubic 
terms in the above expansion are tg and tit2i3- They correspond to Makhhn's 
invariants /i and /12, respectively. 

3. The boundary of the set of separable states 
Using Makhlin's notation, we write the components X, Y, Z in Eq. (|l.l|l as 

X = -^^Sicr^, Y = -^p^a,, Z = ^ /3,j cr^ cr^ , 

where the cr,:'s are the Pauli matrices 

Ol\ f -i \ f I 



In [SJ Table 1] Makhlin lists his 18 invariants /i, . . . , /ig. We shall need only 9 of 
them: 

h = dot/?, h - tr(/3^/3), /3 = tr(/3^/3)2, 

I12 — S/3p, /i4 = eijkeimnSiPlPjmPkn, 

where e^fe is the Levi-Civita symbol. 

In a recent paper jjj Augusiak et al. have shown that S is exactly the set of 
all states p satisfying the inequality detp^ > 0, where T is the operator of partial 
transposition, say with respect to the second party. Denote by J) the real algebraic 
hypersurface in the affine space (l/4)/4 + -ff4,o given by the equation detp = 0. Its 
image under F, which we denote by V^ , is defined by the equation detp^ — 0. 

Let us recall that 5A4 consists of all p > (with trp = 1) having at least one 
zero eigenvalue, and so dA4 C T>. This inclusion is in fact proper, i.e., dA4 is only 
a small portion of the entire hypersurface V. By the Peres-Horodecki criterion |S] 
we know that S ~ Ai O A4^ . Hence, dS is the union of two pieces: 

dS^{Vn M^) U (P^ n M) . 

The first of these pieces has been mentioned in the paper J^l of Vcrstraete et al. 

The polynomial function det p^ is an invariant of G. We have found the following 
expression for it in terms of the above Makhlin's invariants: 

"^""^ ^^ ^ ^' h^^^^ + ^" + ^'^ + ¥^^' + ^^'^ + 1^ ^^^^' " ^^^' 

-16/8 - 16/14 - 16/2^ + lf + I^ + 8/2/4 + 8/2/7 - 2/4/7) . 

There are many papers devoted to the study of the geometry of the sets M and 
S, e.g. ISlElITniiniiniini. it may be of interest to study the boundary dS and 
the determinantal hypersurface T) in more detail. Let us say that a point po € V is 
smooth if the gradient Vdetp does not vanish at po, and otherwise it is a singular 
point of T>. 

The next proposition is valid for any n-dimensional complex Hilbert space Ti.. In 
this more general setting both V and V^ are hypersurfaces in the real affine space 
(l/n)/„ + F„,o- 
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Proposition 3.1. A point po G I) is a singular point of the hypersurface V if and 
only if po has at least two zero eigenvalues. 

Proof. Let us arrange the eigenvalues of po in increasing order 

n-l 

Ai < A2 < • • • < A„ = 1 - ^ Afe 

k=l 

and observe that we must have A„ > 0. Since the unitary group U(n) preserves V 
and maps singular points to singular points, we may assume that po is in fact the 
diagonal matrix 

Po = diag(Ai,A2,.. .,A„_i,l - Ai A„_i). 

For an arbitrary p, the determinant expansion has the form 

detp = P11P22 ■ ■ ■ Pn-l.n-l{l - Pll - • • • - Pn-l,n-l) + P, 

where each term of the polynomial P is at least quadratic in the off diagonal entries 
of p. Hence in evaluating the gradient of det p at the point po the polynomial P 
makes no contribution at all, and it suffices to use just the first term of the above 
expansion, which is written explicitly. We deduce that po is a singular point if and 
only if the following equations hold: 

— (po) = A1A2 • • • Afc • • • A„_i(A„ - Afc) = 0, 1 < fc < 71 - 1 

Opkk 

where the hat means that Afc should be omitted. It is easy to see that these equations 
are satisfied if and only if at least two of the eigenvalues are 0. D 

Note that if n = 2 then it is impossible for both eigenvalues to be 0, which means 
that in this special case the hypersurface V is smooth. Indeed, we have dAi — V, 
I? is a shere, and A4 the corresponding solid 3-dimensional ball. 

The following corollary follows immediately from the proposition. 

Corollary 3.2. A point po G T)^ is a singular point ofD^ if and only if Pq has at 
least two zero eigenvalues. 

Let us return now to the case of two qubits. In that case we shall prove that the 
piece T)^ n A^", where A4° denotes the relative interior of A^, contains no singular 
points of the hypersurface V^ , i.e., for each p e T)^ n A^", the operator p^ has 
exactly one zero eigenvalue. 

Theorem 3.3. The piece of the hypersurface D^ contained in the (relative) interior 
of M. is smooth. 

Proof. Assume that a point po £ V^ n A^° is a singular point of V^ . Then pg is a 
singular point of V and so it has at least two zero eigenvalues. We shall now apply 
an argument of Sanpera et al. JT]. By their Theorem 2, the zero eigenspace of pg 
contains a product vector |e, /), and obviously we have 

(e,/|Po|e,/) = 0. 

But this expression is equivalent to 

(e,r|po|e,r) = 0, 

which is impossible since po G A^", and so po > 0. This contradiction proves the 
theorem. D 
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The following corollary is obvious. 

Corollary 3.4. The piece of the hypersurface V contained in the (relative) interior 
of A4^ is smooth. 
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